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Abstract: The process of exclusive elastic vector meson production in deep inelastic 
scattering is investigated within the dipole model framework supplemented by the small x 
evolution. The dipole-proton amplitude is obtained from the nonlinear Balitsky-Kovchegov 
evolution equation with impact parameter dependence. This dipole amplitude is used to 
compute the differential cross section for exclusive production of J/^', (p, and p vector 
mesons. These numerical calculations are compared with the wide range of experimental 
data from HERA. Good agreement between the experimental data and the calculations is 
found. 
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1. Introduction 

Exclusive diffractive vector meson production in deep inelastic scattering, 7*p — )• Vp^ has 
been extensively studied at high energies at HERA electron(positron)-proton collider. Both 
HI [1~3] and ZEUS [4-7] collaborations performed detailed measurements of this process in 
a wide kinematic range. In this process the final state consists of the elastically scattered 
proton accompanied by a vector meson. It provides important information about the 
electron-proton interaction at small x as well as the size of the proton in impact parameter. 
It can be also used to constrain the form of the vector meson wave functions. The main 
conclusions from these measurements can be summarized as follows (for concise reviews 
of the experimental results see for example [8,9]). The cross section for the elastic vector 
meson production increases with the energy for all values of the (minus) photon virtuality 
. This dependence can be phenomenologically parametrized as ~ where W is the 
energy of the virtual photon-proton system in DIS. The exponent 5 depends strongly on 
the scales involved in the process, namely on the and the mass My of the vector meson. 
The exponent 5 shows a marked increase with both My and scales. For example, for 
p production the variation of 5 with is from 0.2 for photoproduction up to ~ 0.8 for 
values of ~ 35 GeV^. This is consistent with the picture of the soft-hard transition 
from small to large scales. The value of the 5 in the photoproduction region, and at low 
(5^, is consistent with the soft exchange, while in the high region perturbative gluon 
exchange is the dominating process. For low values of multiple interactions between the 
quark-antiquark system and a proton also lead to smaller values of 5, see for example [10] . 
For J/^' production the variation of 5 with is negligible in the measured range (from 
photoproduction to electroproduction with ~ 30 GeV^), and its value is about ~ 0.8, 
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indicating that the heavy mass of the J/^' vector meson provides the hard scale in this 
process. 

The measurement of the momentum transfer t in this process provides the information 
about the transverse size of the interaction region, see for example [11]. These measure- 
ments indicate a strong dependence of the differential cross section on t, which can be 
parametrized by an exponential exp{—B£)\t\). The experimental data show that the value 
of Bd decreases from about ~ 10 — 12 GeV~^ at lowest values of {Q^ + My) to about 
5 GeV~^ at scale + My ~ 10 GeV^, and stays constant for larger scales. This can be 
interpreted as follows: the value of parameter Bjj is closely related to the transverse size 
of the interaction region which is a combination (or more precisely, a convolution) of the 
sizes set by the vector meson and the proton. At low values of + My it is the first size 
that prevails, and as a result the value of B^ strongly varies with this scale. At larger 
values of + My the data indicate that the vector meson size is much smaller than the 
typical size of the proton, (or more precisely the radius of the gluon density in the proton), 
and it is this 'gluonic' size of the proton which is dominating Bd and thus visible at larger 
values of + My. This value is about 5 GeV^^ or equivalently ~ 0.6 fm, and it is smaller 
than the electromagnetic size of the proton. The data also indicate a slow variation of the 
slope Bd with the increasing energy of the interaction W. The slope can be parametrized 
by the following Regge inspired formula Bd = Bq + Aa'jpln{W/Wo), with the fitted value 
of a'jp ~ 0.115 ± 0.018(+0.008) (-0.015) [4] and a'jp ~ 0.164 ± 0.028 ± 0.030 [2] for J/^ 
photoproduction (parameter Wq = 90 GeV). This shows that the size of the gluon distri- 
bution in impact parameter space is increasing with rising energy, which is indicative of 
Gribov-type diffusion of gluons in transverse space. 

In this paper we analyze elastic diffractive vector meson production using the dipole 
model approach for small x physics. The novel aspect of our analysis is that the dipole 
scattering amplitude is obtained from the numerical solution to the impact parameter de- 
pendent Balitsky-Kovchegov (BK) equation [12-14]. Numerous studies, which utilized the 
dipole model, were performed up to date whose goal was a description of the exclusive 
diffractive vector meson production at small x, see for example [10,15-20]. As mentioned 
above, the t-dependence of the vector meson production cross section is particularly in- 
teresting as it is directly related to the transverse size of the interaction region. In the 
framework of the dipole model this is encoded in the dipole-proton scattering amplitude 
and its dependence on the impact parameter. The increase of the t slope, Bd, with the 
energy W should be directly coupled with the increase of the transverse interaction area 
caused by the diffusion of the dipoles in coordinate space. This diffusion in transverse 
coordinate space is present in the dipole branching Monte-Carlo [20-23] as well as the 
Balitsky-Kovchegov equation when the impact parameter dependence is taken into ac- 
count [24,25]. The main goal of the present work is to use the numerical solution to the 
BK equation, which includes the full impact parameter dependence, and verify its com- 
patibility with the experimental data on the exclusive vector meson production in deep 
inelastic lepton-proton scattering. 

The outline of the paper is the following, in the next section we recall the formalism of 
exclusive vector meson production in the dipole model, in Subsec. 2.1 we discuss the basic 
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Figure 1: Schematic representation of the inclusive cross section for deep inelastic scattering in 
the dipole model. Virtual photon fluctuates into qq pair and interacts with the hadronic target. 
The dynamics of the interaction is encoded in the dipole - target scattering amplitude N. 



facts which pertain to the solution of the BK equation with impact parameter dependence, 
and discuss modifications due to confinement. In Subsec. 2.2 additional modifications 
are outlined, such as the non-perturbative modification of the photon wave function and 
inclusion of the skewed effect into the initial gluon distribution, as well as inclusion of the 
real part of the scattering amplitude. This is followed by the description of the model for 
the vector meson wave function which is used for the calculation. In Sec. ^ we present our 
results, which include the comparison of the theoretical calculation with a wide range of 
experimental data from HERA on p, <j) and J/^ production. We compute the cross section 
integrated over the momentum transfer t and investigate its W and dependence. Then 
we compute the differential cross section with t dependence in bins of W and Q^. The 
experimentally measured scale and energy dependence of the slope Bd is also compared 
with theoretical calculations. Finally, in the last section we state the conclusions and 
present an outlook. 



2. Exclusive vector meson production in the dipole model 

The dipole model [26, 27] is a very useful tool in evaluating many processes at small values 
of X. One of the advantages of this approach is the possibility of including the multiple 
parton scattering effects. It has been originally formulated for the description of deep 
inelastic lepton-proton (or nucleus) scattering at small x. In this picture, utilizing the 
leading logarithmic approximation in x, the incoming electron emits a virtual photon which 
fluctuates into a quark-antiquark pair (a dipole). This color dipole then subsequently 
interacts with the parton constituents of the nucleon, as is illustrated in Fig. |l[ 

The interaction of the dipole pair with the target is given by the scattering amplitude 
A^. This quantity is non-perturbative in principle, however its small x evolution can be 
found from the BK equation. The qq pair is characterized by a dipole size which is defined 
as a separation distance of the color charges xqi = xq — xi (where xq and xi are the 
positions of the q and q in transverse space). ^ The transverse momentum of the quarks in 

""^In this paper we shall denote vector quantities in bold, otherwise they should be read as magnitudes 
of the associated vector. Also, alternatively we will be also using here the notation for the dipole size to be 
r = a;oi and impact parameter b — isiill. 
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the dipole is on the order of ~ — where large dipoles correspond to the infra-red region 
and need to be regulated. The scattering amplitude A^(r, b; Y) contains all the information 
about the dynamics of the strong interaction. It depends on the dipole size r, on the impact 
parameter of the dipole with respect to the target b and the rapidity Y. In the following 
analysis the full dependence of the scattering amplitude on the impact parameter b will 
be taken into account. 

The structure functions F2 and Fl for the proton can be evaluated using the following 
standard formulae in the dipole picture in the transverse coordinate representation 



47r^Q;em J Jo 

and ^ ^ 

Fl{Q^x) = — ^ / d^r [ dz|^L(r,z,Q2)|2 (r,x) . (2.2) 

The dipole-proton cross section cjdip can be obtained from the scattering amplitude by 
integrating over the impact parameter b 

(Tdipir, ^) = 2 y ^(r, b; F) , Y = ln 1/x . (2.3) 

Since the amplitude is dimensionless and the integration is over the impact parameter, 
the dipole cross section fidip has obviously a dimension of the area. We see therefore that 
although the inclusive quantities are sensitive to the size of the interaction area, the details 
of the impact parameter profile are not directly accessible through this process. 

The quantities '^{r,Q'^ ,Y)x/l are the photon wave functions. They describe the dis- 
sociation of a photon into a qq pair and can be calculated from perturbation theory. The 
photon wave function has the following form for the case of transverse photon polarization 

l^Tir, z, Q')\' = ^ E 4 + (1 - ^f] Q}Ki (Qf^) + {Qfr)) , (2-4) 



and for longitudinal polarization 



^Ur,z,Q'r = ^Y.^}i^Q"'"(^-')"^o{Qfr)) ■ (2.5) 



In the above equations = z{l — z)Q'^ + m^, where —Q^ is the photon virtuality and 
z, (1 — z) are the fractions of the longitudinal momentum of the photon carried by the 
quarks. In addition Kq^i are modified Bessel functions of the second kind. The summations 
are over the active quark flavors /, of charge ej, and mass rrif. 

The dipole picture can be also used to compute diffractive processes. Here, we are 
interested in the process of the exclusive, diffractive production of the vector-meson 7*p — )• 
Vp' . The amplitude for this process is schematically illustrated in Fig. |2|. The virtual 
photon still fluctuates into qq pair, which then interacts with the proton and a vector 
meson is formed, which is measured in the final state. The proton scatters elastically. 
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Figure 2: Schematic representation of exclusive vector meson production in the dipole model at 
small X. A virtual photon fluctuates into a qq pair and interacts with the target (proton). After the 
interaction the vector meson V is formed which is measured in the final state. The proton scatters 
elastically with some momentum transfer t. 



its 4- momentum in the initial state is p and in the final state is p' . The formula for the 
amplitude for this process reads 

A{x,A,Q) = J] /c?'r / dzvI/ft,;,*(r,z,Q2)AA(a,,r,A)vI'^^^.(r,z) , (2.6) 
h.h •' •' 

where h{h) is the helicity of quark (antiquark) and ^^^^ is the vector meson wave function. 
A is the 2-dimensional momentum transfer related to the Mandelstam variable t = — A^. 
The differential cross section for the process is given by 

f = ^|A(.,A.Q)r (2.7) 

The amplitude M{x, r, A) can be related to the scattering amplitude N{x, r, b) introduced 
earlier, the amplitude in the impact parameter representation through the appropriate 
2-dimensional Fourier transform 

M{x, r, A) = 2 j (fh N{x, r, b) e'^'^ . (2.8) 
In this notation the dipole cross section, (compare (|2.3|) ), is 

adipix, r) = Im iM{x, r, A = 0) , (2.9) 

which is the expression for the optical theorem for scattering of dipoles. 

This process, through its dependence on the momentum transfer t, offers a unique 
possibility of constraining the impact parameter profile of the dipole scattering amplitude. 

Formulae (^^) and ( |2.8D were original expressions derived under the assumption that 
the dipole size is much smaller than the proton. In Ref. [28], a correction due to the finite 
size of the dipole was calculated. It was shown that in the non-forward case, A 7^ 0, the 
amplitude can be written in the similar form as above with the modification of the ( p.^ ) 
to include the exponential factor exp(— i(l — z)r ■ A) in the following way 

Af{x, r, A, z) = 2 J d^h N{x, r, b) e^^-(b-{i-^)r) . (2.10) 
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This modification was included in the calculation [29] and it was shown that it has a non- 
negligible effect on cross sections, especially on the values of the Bd slope which controls 
the t-dependence as a function of the scale + My . 

2.1 Dipole scattering amplitude from impact parameter dependent BK evolu- 
tion 

The dipole-proton scattering amplitude iV(r, b; Y) at high values of rapidity Y (or small 
x) is found from the solution to the nonlinear integro-differential Balitsky-Kovchegov (BK) 
evolution equation [12-14,30]. The BK evolution equation can be represented in the fol- 
lowing form: 



dY 



/d X2 
-^/C(xoi,Xi2,2;o2;as,m) [iVxoxa + ^xaxi -^'xoxi - iVxoX2^^X2Xi] • (2.11) 



In the above equation we used the shorthand notation for the arguments of the ampli- 
tude iVxjXj = N{rij = Xj — Xj,bjj = i(xj -|- Xj);y) which depends on the two transverse 
positions Xj and Xj and on the rapidity Y. The branching kernel IC{xoi, xu, xq2; Osjin) 
depends on the dipole sizes involved and contains all information about the splitting of the 
dipoles. In addition, it depends on the running coupling a^. The way the strong coupling 
runs will be specified later in this work. We have also indicated that the kernel depends 
on the infra-red cutoff m which we impose in order to regulate large dipoles. 

Eq. ( p. 11 ) is a differential equation in rapidity and hence suitable initial conditions need 



to be specified at some initial value of rapidity Y = Yq. As in the previous work [31] we are 
choosing to use the initial condition in the form of the Glauber - Mueller parametrization 
with (most of) the parameters equivalent to those used in Ref. [29] 

Ncuir, b;Y = In 1/x) = 1 - exp (^-^r^xg{x, 7f)T{h)^ , (2.12) 

with ^ 

T{h) = i-e^ . (2.13) 



In formula ( p. 12 ) the function xg{x,'rf') is the integrated gluon density function and 



T{h) is the density profile of the target in transverse space with the extension set by the 



parameter Bg- The integrated gluon density in (2.12) was also taken from fits performed 
in [10]. Scale parameter in the gluon density is set to he = ^'^ + ^ with parameters /io 
and C = 2 set to obtain the best description of the data. The values of these parameters 
are given in Table |^. We use ( 2. 12] ) as the initial condition at Iq = IhI/xq, xq = 



and evolve the amplitude with the BK equation to obtain the solution at lower values of 
x < xq. We also note that the initial condition ( 2.12| ) depends only on the absolute values 



of the dipole size and impact parameter. A nontrivial dependence on the angle between 
vectors r and b is not present in the initial condition, instead being dynamically generated 
when the initial condition is evolved with the BK equation. 

The BK equation was solved numerically by discretizing the scattering amplitude in 
terms of variables (logj^g lo§io where 9 is the angle between the impact parameter 
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b and the dipole size r. The amphtude N{r, b, cos 0) was placed on a grid with dimensions 
200r X 200ft X 20g. More details about the techniques and the properties of this solution 
can be found in Refs. [24,25]. 

The running of the coupling in the BK kernel is a next-to- leading effect which has been 
evaluated in [32,33]. In this calculation we utilize the prescription from Ref. [32] which is 
of the form 



/C = as(xoi) 



1 



02) 



'-02 \'^s(^12) 

In the above equation, we use 



1 + 



1 



'12 



I2) 



a 



^(^02) 



1 + 



12-^02 



(2.14) 



a. 



h In 



A, 



QCD 



(2.15) 



'^'^12tt^ ^'^^ ?T.f is the number of active flavors. The /x parameter effectively 



61n[A-V^l 



In our simulations we used 
0.246 GeV, and 



Here h = — y2tt^ '^'^^^ 
freezes the coupling at large dipole sizes at a^^fr 

/X = 0.52 GeV as an infra-red regulator for the strong coupling, Aqcd 
the rescaled strong coupling is defined as as = ""^^ ■ 

The evolution equation ( |2.11| ) has been derived in perturbation theory and it does not 
include any effects of confinement. The branching kernel is power-like in the dipole sizes 
and it allows for the splitting of a parent dipole into pair of arbitrarily large daughter 
dipoles. As a result the impact parameter dependence of the amplitude contains Coulomb- 
like power tails in impact parameter, which need to be regulated. This has to be done by 
cutting off the large dipole sizes in the branching kernel. This is done by including a mass 
parameter, m, which accounts for the effect of color confinement. In previous works we 
have tested several prescriptions and found the simulation with the simple cutoff using the 
theta functions gave best description of the inclusive data. We have therefore used this 
scenario for the calculation in this paper. The modified kernel which has been used for the 
calculation of a solution has therefore the form of the kernel ( p.l4| ) with theta functions 



1 



-02 



as{xl2) 

as{x'^ 



12) 



1 + 



x 



12 



as{xl2) 

as{x'^ 



02) 



1 + 



^12^02 



X 0(^:2 



^02 



^12) 



(2.16) 



where m is the mass regulator. The value of the parameter m has been fitted to obtain 
the best description of the data. Note that the similar procedure has been utilized in 



Refs. [20,22], with the value of the parameter 



to be equal around 3 GeV . The 



formula (|2.12| ) is used as an initial condition for the BK evolution for dipoles with sizes 
smaller than the cutoff 1/m, for dipoles larger than the cutoff the initial condition is set 
to zero. 



2.2 Corrections to the dipole scattering amplitude and photon wave function 

There are several additional phenomenological corrections which we have included in this 
calculation. As we will see explicitly they have a non- negligible impact on the calculations. 
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The first of these is to take into account the effect of the skewed gluon distribution. This 
is necessary as the two gluons exchanged in the production of the vector mesons need not 
have the same longitudinal momentum fractions x and x' . The skewed effect vanishes at 
small X (in the leading logarithmic limit in In 1/x), however it can be substantial correction 
when the energy is not very large. As suggested in Ref. [34] this correction can be taken 
into account by a multiplicative factor on the standard gluon distribution as follows 

2,, , 2(^^-^-+3)r(A,fc + 5/2) 

{xg{x,r] )),k = xg{x,r] ) ^ — — — — , (2.17) 

Vvr 1 [Ask + 4j 

where 

_ dln{xg{x,r]'^)) 

= 51n(l/x) • (^-^^^ 

Strictly speaking, in the dipole formalism there is no integrated gluon distribution, but 
rather the dipole scattering amplitude. However, since the initial condition is taken as the 



Glauber-Mueller model in the form (2.12), with explicit dependence on the gluon density 
xg{x, "rf), we can implement the correction for skewedness inside the initial condition. Then 



the evolution to lower values of x is computed according to Eq. ( 2.11| ). 



In addition we have included the correction for the real part of the scattering amplitude. 
This effect can be taken into account by multiplying the amplitude by (1 + /3^) where 

(3 = tan(7rA^/2) , (2.19) 

with 

8MX<'V^'') 

^ = SHIM ■ <^-^»' 

In the above equation (3 is the ratio of the real to imaginary part of the scattering amplitude. 
A similar procedure was used previously in the other descriptions of the vector meson 
production [17,19,34]. 

Finally, a correction to the photon wave function is necessary in order to modify and 
enhance the contribution at low Q^. The photon light-cone wave function given in Eqs. 



( |2.4| ) and (2^) has been derived in the perturbation theory under the assumption of the 
presence of hard scale , which in turn leads to the small dipole sizes r. This statement is 
most accurate for the longitudinal contribution, but in the case of the transverse contribu- 
tion the endpoint singularities in z cause the distribution in dipole sizes to be broader, even 
for very large scales Q^. At low values of photon virtuality the corresponding dipoles 
which contribute to the cross section are very large and non-perturbative. Consequently 
the cross section with which they interact with the proton is also large. Therefore, at these 
large dipole sizes the photon has a hadronic component, which is non-perturbative. To ac- 
count for this effect we use a modification of the photon wave function in a way suggested 
in Ref. [19] and also used in Ref. [20] 



\^^\^ — 2752 • (2-21 
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1.16547 GeV 




Bg 


3.65 GeV-2 


D 


2 




m 


0.37 GeV 




0.52 GeV 




R 


6.8 GeV-i 




2.55042 




B 


6.0 




0.01980 




u? 


0.2 GeV^ 



Table 1: Values of the free parameters used in the calculations. 

The constants i?, uj and R are the parameters which have to be adjusted to fit the 
data. This factor provides an enhancement for dipoles which have a hadronic size. The 
numerical values of the parameters B,uj,R which are used in the calculation are given in 
Table |l[ 

2.3 Vector meson wave function 

Several different models for the vector meson wave function exist in the literature, for ex- 
ample see [15-17,35-38]. Typically, there are many uncertainties in obtaining the wave 
functions for the vector mesons, however they are constrained by model independent fea- 
tures. First of all, il^y^ has to satisfy the following normalization condition 

1 = E / d''rdz\ij'^\z,r)\'' . (2.22) 
h'h 

In addition, the value of the wave function at the origin is related to the leptonic decay 
width T{V — )• e~^e~) of the vector meson, given by the following formula 

^'dzVy(^,r = 0) = y^^, where (0| Je^^(0)| V) = e,/y my e'^ . (2.23) 

Here fv is the coupling of the meson to the electromagnetic current and ey the isospin 
factor, which is the effective charge of the quarks in units of the elementary charge e: for 
the p meson it is the charge of the combination {uu—dd)/y/2, i.e. ey = l/\/2. Finally, one 
requires that the mean radius be consistent with the electromagnetic radius of the vector 
meson. 

In this paper we utilize the model proposed in [15, 17]. In this approach the information 
from spectroscopic models is used to constrain the long distance physics. One assumes that 
the meson is composed of a constituent quark and antiquark which move in a harmonic 
oscillator potential. This results in a wave function which has a gaussian dependence on 
the spatial separation between the quarks. Additionally, this model is supplemented by 
the short-distance physics driven by QCD exchange of hard gluons between the valence 
quarks of the vector meson. Finally, a relativization technique has to be applied to the 
wave function. 

The wave function of the vector meson is given by 

i^f\z,r) = V-fe _ - ^) - V.^ + ml]<P{r, z) . (2.24) 
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P 


</> 




*o(15)t 


0.025 


0.028 


0.039 


^o(15)l 


0.024 


0.025 


0.039 




12.77 


11.0 


2.188 



Table 2: Values of the free parameters used in the vector meson wave functions. 



In this model, the radial wave function 0(r, z) appearing in eg. ( p. 24 ) satisfies the 
following normalization condition: 

1 = ?^ r^T;^ [d'rlmyiz,r) + iz' + il-zf)idrHz,r)A . (2.25) 



2tt Jo z^il - zY 



This equation is essentially the normalization condition ( p. 22 ) applied to the transversely 



polarized meson wave function. The function </)(r, z)(/,/j^) is defined by: 



r , / miW \ ( 2z{\-z)r^\ / m^iJM 1 

X |4.(1 - .)^/^exp y-^) e>=P (-^Ip^) »P (^j } ^ P-^") 

The masses of the light quarks are taken to be 0.14 GeV and the charm mass is taken to 
be 1.4 GeV. 

The other two parameters ^'ol/t(1'S') and i?^ are chosen by taking into account several 
constraints: the normalization condition ( |2.25| ) for the wave function has to be satisfied, 
and the value of the leptonic decay width must agree with the experimental measurement. 
Additionally, the mean radius of the meson has to be of the order of a hadronic scale. The 
normalization 'I'ol/t(15') is different for the longitudinal and transverse wave functions of 
each vector meson. On the other hand the value of -R^ is the same for each polarization but 
different for each vector meson. The values of the parameters used for the vector meson 
wave functions can be found in Table |2[ 



3. Comparison with the experimental data 

In this section we compare the results of our calculation based on the numerical solution of 
the BK equation with the impact parameter dependence to the data from HI and ZEUS 
on /), 0, and J/^ production. 

As discussed above, the initial condition for the evolution in x is given by Eq. ( 2.12| ). 



This formula contains several free parameters, such as the value of Bq, the parameters 
Ag^\g in the initial gluon distribution xg{x,rf') = AgX^'^o{l — x)^'^, and the parameters 



//Q and D which enter into the definition of the scale rf' = + /ig- The initial impact pa- 
rameter profile T{h) (|2.13| ) is assumed to be gaussian. This profile is substantially changed 
during the course of the evolution with x. In addition, the kernel contains two scales: 
the scale in the running of the coupling as well as the mass parameter m. All these 
parameters enter into the calculation of the dipole scattering amplitude which is found by 
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the iterative solution of the Eq. ( |2.11| ). The computational time that is required to obtain 
the solution for each set of parameter values is rather long (on the order of 24 hrs on 32 
cores). Therefore we have chosen to vary only a subset of these parameters. Parameters 
A, \g,C, Ho were fixed and their values were taken from Ref. [10]. The only parameter 
which was varied in the initial condition Eq. ( 2.12| ) was Bq- The value of Bq was adjusted 
to obtain the best description of the t dependence of the differential cross section da/dt. In 
particular this parameter controls the magnitude of the experimentally determined value 
Bd, which is the slope of the t dependence, i.e. da/dt ~ exp(— In addition, the 
mass parameter in the kernel m, which controls the W dependence of the slope B^ was 
also fitted, together with Bg in Eq. ( |2.12 ), to obtain the best description of the data. In 
our analysis we set m = ^^jg— but these two parameters need not be directly linked. 

There are also parameters that are present in the formulae for the cross section for 
the production of the various vector mesons. These include the quark masses as well as 
the parameters in the vector meson wave function and non-perturbative correction to the 
photon wave function. The variation of these parameters does not cost so much computa- 
tional time therefore a larger number of them were varied for the best fit. The complete 
set of parameters used in this paper can be found in Table ||. 

Using the framework described above we have performed a fit of the calculation to 
the inclusive data on the structure functions in DIS. Note that this was necessary as the 
details of the procedure presented above are slightly different than that originally used 
in [31]. Good agreement between the inclusive data and the calculation was found using 
this procedure. 



3.1 Cross sections for exclusive vector meson production integrated over t 

Let us first show the comparison between the calculation based on the dipole model with 
BK equation and the experimental data on the cross section for the process of exclusive 
diffractive electroproduction of vector mesons, where the cross section has been integrated 
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Figure 4: Cross section (t{Q^,W^) for tlie vector meson production plotted as a function of 
{Q^ + My) for (j) (left plot) and J/ijj (right plot) elastic production. The experimental data are 
from [2,3,5,6]. 



over the momentum transfer t. The experimental data from HI and ZEUS for p [1,3,7], 
(f> [3,6] and J/^' [2,5] were used. 

Figures 1^ and ^ shows the cross section for production of p, cj) and J/\I' vector mesons as 
a function of variable (Q^+My), where My is the mass squared of the corresponding vector 
meson. This variable is commonly used instead of itself as it provides the scale for the 
vector meson. One cannot, however, expect that cross sections for different vector mesons 
will behave identically when this variable is used. We therefore show the cross sections for 
different species of vector mesons separately. The momentum transfer t has been integrated 
in the ranges provided by the experimental data. For the data on p production from [1] 
the range in t is \t\ < 0.5 GeV^ and \t\ < 1.0 GeV^ for data from [7]; \t\ < 3 GeV^ for 
data from [3]. For the data on (j) production from [6] the range is |t| < 0.6 GeV^ and 
|t| < 3 GeV^ for data from [3]. Finally for J/"^ the range of \t\ < 1.2 GeV^ in [2] and 
jtj < 1.0 GeV^ for data from [5]. The theoretical curves shown in Figs. |3|,§ are evaluated 
for constant energy W where W = 90 GeV for the J/^ cross section and W = 75 GeV or 
W = 90 GeV for the (j) and p cross sections depending on the data set used. 

We have tested the sensitivity of the cross section due to the inclusion of the various 
corrections discussed in Sec. 2^. This is shown in two plots in Fig. |5|. In the left plot 
the effect of including the correction due to the skewed effect in the initial condition is 
shown with data for J/^ production. We see that the effect is the substantial change 
in the normalization of the cross section, without distortion of the dependence, and 
calculation without skewed effect correction completely fails to match the data. We note 
that even with the substantial variation of the free parameters it would be difficult to match 
the normalization of the cross section without this additional effect. 

In the plot on the right hand side we demonstrate the effect of including the non- 
perturbative modification of the photon wave function in the form of Eq. ( 2.21| ). This 
correction is especially important for low values of and light vector mesons. This is to 
be expected as this correction vanishes for large in order to reproduce the perturbative 
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Figure 5: Illustration of the effect of various corrections on the integrated cross section for the 
vector meson production. Left: the solid line represents the calculation without skewedness and the 
real part of the scattering amplitude, while the dashed line represents calculation which includes 
skewedness and a correction to account for the real part of the scattering amplitude in the initial 
condition. No correction to the photon wave function was made on either of these curves. Right: 
both curves have the initial condition corrected for the real part of the scattering amplitude and 
include a skewed gluon distribution correction. The solid line represents the inclusion of the photon 
wave function correction and the dashed line is the calculation without the photon wave function 
correction term. 



expression for the photon wave function. Thus this correction has a negligible effect at 
large values of > 20 GeV^ and the cross section for J/^* production. 

The W energy dependence of the cross sections is shown in Figs. ^,|^J8| for /),(/), J/*? 
respectively. The different curves are plotted in different bins of Q^. Overall trend in both 
cases, i.e. for and W dependence, is such that the calculations describe the dependence 
in and W dependence very well for the case of (p and J/^'. The data for p production 
are not well described, in particular the normalization in this case is systematically little 
bit low, especially for lowest values of Q^. This region is however the one which is not 
under perturbative control and some unknown non-perturbative corrections, related for 
example to the exact form of the wave function, may play an important role in this region. 

The ratio of the longitudinal to the transverse part of the cross section R = ^ was 
analyzed as well. This ratio has a significant sensitivity to the exact form of the wave 
functions used. In Fig. ^ the calculation is compared with the experimental data as a 
function of Q^. The data sets shown in figures are for very wide bins on W, and therefore 
we have shown the curves which correspond to the middle value of the bins. 

By inspecting the formulae for the transverse and longitudinal cross sections one would 
think that the ratio should be approximately independent of the energy W. This is not 
entirely true as the longitudinal and transverse components of the cross sections are sen- 
sitive to somewhat different distributions of the dipole size configurations in the photon 
wave functions. It is well known that the transverse photon wave function takes more 
contributions from the large dipole sizes, and therefore one should expect that the energy 
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Figure 6: W dependence of the vector meson cross section for elastic production of p. The 
experimental data are from [1,3,6,7]. 

fp=>*p t'p=>*p 




W (GeV) W (GeV) 

Figure 7: W dependence of the vector meson cross section for elastic production of (p. The 
experimental data are from [3,6]. 

dependence of or should be flatter than that of ai. Experimental data show however that 
the ratio R does not depend on the W indicating that the large dipole components may 
be suppressed in the transversly polarized exclusive vector meson production. The overall 
description is very good with the exception of the (j) data where the theoretical curves 
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Figure 8: W dependence of the vector meson cross section for elastic production of J/'^. The 
experimental data are from [2,5]. 



Y' P => P P T' P -> P P 
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Figure 9: The ratio of i? = ^ cross section for longitudinally to transverly polarized vector 
mesons. The plots are for p (top left W = 75 GeV and right = 90 GeV), (j) (bottom left 
W ^90 GeV), and J/^- (bottom right W ^ 75 GeV). 



indicate the growth with the energy which is slightly too fast as compared with the data. 
This trend is similar to what was observed in [10] which indicates a generic feature of the 
calculations based on the dipole model. 
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3.2 Differential cross section and t-distribution in vector meson production 

The t-distribution of the differential cross section for elastic vector meson production is 
key in unraveling the impact parameter profile of the target at small-x. The t-dependence 
can be related to the impact parameter dependence via two-dimensional Fourier transform 
of the amplitude as indicated in Eqs. ( |2.6| ),(^I^). The differential cross section is usually 
parameterized as ^ oc e"^^'*' in bins of and W . The dimensionful slope parameter By) 
thus contains the information on the spatial distribution of the interaction region in the 
scattering process. Three plots in Figure |l^ show the dependence of the slope parameter 
on the variable + My for 0, and Jj'^. The theoretical curves follow the trend of 
the experimental data. We observe that for the p production the dependence of By) on 
is well described but the normalization is underestimated, which is most probably 
related to the lower normalization for the resulting integrated cross section. The decrease 
of the slope for low values of + My is related to the initial dependence on the size 
of the vector meson. For larger values of the dependence flattens to a common value 
of By, ~ 4 GeV^^ for each of the vector meson species. This flattening and universality 
at large values of indicates that in this regime the By) indeed characterizes the size 
of the proton through the interaction with the small probe which is high dipole. This 
characteristic size of the gluon density inside the proton \J (r^) ~ 0.6 fm is markedly smaller 
than the electromagnetic radius which is of the order ~ 0.8 fm. This indicates that the 
gluon distribution differs from the spatial extension of the quarks in the proton. 

In Fig. |ll| we show the same quantity Bd but as a function of W for two different 
values of for both J/'f and p. While the error bars on the experimental data for Bd 
are relatively large, we see that the theoretical curves describe very well the increasing 
trend of the data, which is especially visible in the bin for lower value of (in fact the 
bin with higher is consistent with the flat dependence as well). In the calculation 
presented, the energy dependence of the slope is naturally obtained from the evolution of 
the dipole scattering amplitude with the energy. The diffusion of the dipoles in the impact 
parameter space is what provides the change of B^ with energy. Since it is encoded in 
the BK evolution it naturally leads to the broadening of the impact parameter profile with 
the energy. The normalization and slope (in energy) of the B^, however, do depend on 
the two free parameters in the calculation which are not calculated from first principles. 
The intercept depends on the value of Bq which is used in the initial condition provided 
by Glauber-Mueller formula, Eqs. ( |2.12| ), ( |2.13| ). The slope in energy of Bd is, on the 
other hand, controlled by the value of the mass m in the dipole evolution kernel which 
cuts off the large dipole sizes. The Bd parameter for J/^ is very well described both in 
normalization and in the W dependence. On the other hand the the normalization of p 
is once again low, nevertheless the dependence in W is well described. We note that the 
value of Bg was fitted from the normalization of the Bd slope of J/^ and that the mass m 
parameter which is related to the W dependence is correlated in the calculation with Bq. 
The comparison of the calculation with the data shows that perhaps the initial size in the 
dipole scattering amplitude may not be universal between p and J/^ production or that 
the wave function of p needs to include additional non-perturbative components which will 
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Figure 10: Dependence of the slope parameter Bd on combined variable + My for p, 0, J/^ 



increase the interaction size for this meson. 

Finally we note that, the differential cross section of J/^ production compares very 
favorably with the HI data [2] both in t and in W dependance, which is shown in Figs. (|l2|) 
and (pi). 



4. Conclusions 



In this paper we have computed the integrated and differential cross section for exclusive 
diffractive vector meson production in deep inelastic scattering using the dipole model 
framework for small x. The dipole - target scattering amplitude was obtained from the 
numerical solution to the impact parameter dependent Balitsky-Kovchegov equation. We 
have found that the overall description of the data is very good, meaning that the calcula- 
tion based on BK equation with impact parameter dependence is able to reproduce all of 
the features seen in the data. Several comments are in order: 

1. For the good description of the experimental data within the framework presented it 
is important to include additional corrections. Non-perturbative modification of the 
photon wave function at low values of was necessary, which enhanced the cross 
sections of the dipoles with sizes of the order of the hadronic scales. 
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Figure 11: Dependence of the slope parameter Bo versus W for J /^jj and p production. Data are 
from HI experiment [2] and [7]. 

2. The skewedness effect was included in the gluon density distribution. This distri- 
bution is present in the initial conditions for the small x evolution. This had a 
substantial impact on the normalization of the resulting cross section and helped to 
bring the calculations to agreement with the experimental data. 

3. The BK equation was modified to include confinement effects by cutting off large 
dipole sizes. The parameter rmax = which sets the maximal size of the interac- 
tion, together with the initial proton size control the slope of the differential cross 
section with respect to t as well as its variation with the energy. The presented cal- 
culation shows very good agreement with the experimental data on B^, including its 
W dependence in the case of J/^. The slope of Bd is reproduced for p but the nor- 
malization remains low. The W dependence is generated dynamically in the dipole 
evolution. The speed of this increase is controlled by the parameter rmax = ^ which 
is not calculable from perturbation theory and needs to be adjusted. 

4. The calculation presented includes the running coupling in the evolution, but misses 



-18- 




Figure 12: The differential cross section of J/^ production as a function of W for fixed in bins 
of momentum transfer t, data from HI [2]. 




Figure 13: Differential cross section of J/^E* production for a fixed W in bins of as a function 
of momentum transfer \t\. Calculations were done with W = lOOGeV and W — 90GeV. The 
experimental data are from HI experiment [2]. 



other important NLL effects which are known to be non-negligible. These should 
help to bring the calculation to a better agreement with the data, especially as far 
as the W dependence is concerned. The analysis which includes these effects is thus 
left for further investigation. 
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